Astronomy & Astrophysics manuscript no. Paper 


5th Februai-y 2008 


(DOI: will be inserted by hand later) 





Variability and periodicity of field M dwarfs revealed by 

multichannel monitoring 

Boris Rockenfeller, Coryn A. L. Bailer- Jones*, and Reinhard Mundt 

Max-Planck-Institut fiir Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 
Received 6 September 2005 / Accepted 15 November 2005 

Abstract We present simultaneous, multiband photometric monitoring of 19 field dwarfs covering most of the M spectral 
sequence (M2-M9). Significant variability was found in seven objects in at least one out of the three channels I, R and G. 
Periodic variability was tested with a CLEAN power spectral analysis. Two objects, LHS370 (M5V) and 2M1707+64 (M9V), 
show periods of 5.9 ± 2.0 and 3.65 ±0.1 hours respectively. On account of the agreement with the typical values of vsini 
published for M dwarfs (Mohanty & Basri 2003), we claim these to be the objects' rotation periods. Three further objects show 
possible p eriods of a few hou rs. Comparing the variability amplitude in each channel with predictions based on the synthetic 
spectra of lAllard et alJ i200ll ). we investigated the source of variability in LHS370 and 2M1707+64. For the latter, we find 
evidence for the presence of magnetically-induced cool spots at a temperature contrast of 4 - 8 %, with a projected surface 
coverage factor of less than 0.075. Moreover, we can rule out dust clouds (as represented by the COND or DUSTY models) 
as the cause of the variability. No conclusion can be drawn in the case of LHS370. Comparing the frequency of occurrence 
of variability in this and various L dwarf samples published over the past few years, we find that variability is more common 
in field L dwarfs than in field M dwarfs (for amplitudes larger than 0.005 mag on timescales of 0.5 to 20 hours). Using the 
homogeneous data sets of this work and Bailer- Jones & Mundl (2001), we find fractions of variable objects of 0.21 ±0.11 
among field M dwarfs and 0.70 ± 0.26 among field L dwarfs (and 0.29 ± 0.13, 0.48 ± 0.12 respectively if we take into account 
a larger yet more inhomogeneous sample). This is marginally significant (la deviation) and implies a change in the physical 
nature and/or extent of surface features when moving from M to L dwarfs. 
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1. Introduction 

Photometric and spectroscopic monitoring of very low mass 
stars and substellar objects has been performed over the last 
ten years by several groups and has led to new insights into 
the processes occurring in these objects. Since stars on the 
main sequence become fully convective below a mass of about 
0.3 - 0.4 Mq (roughly coiTesponding to a spectral type of M3 
to M4), magnetic fields cannot be maintained in stars of lower 
mass or later type via an o-Q dynamo. An alternative dynamo 
may come into operation, for example an a^-dynamo (Chabrier 
& Kiiker astro-ph/05 10075). A change in the dynamo mech- 
anism with spectral type could directly affect the activity of 
M dwarfs as measured by Hq, emission, for example. Early 
M dwarfs in most cases do not show any detectable Hq, emis- 
sion whereas it increases at mid and late M type. Furthermore, 
a saturation-type relation between rotatio n a nd activity has 
been c onfirmed by Mohantv & Basril (l2003l) and lPelfosse et alJ 
( Il998 l). The detection of flares on these objects also proves 
the existence of magnetic fields. However, magnetic activity 
strongly decreases again with early L type dWest et al.ll2004l) 
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which can be explained via an incr easingly neutra l photosphere 
and the onset of dust formation ( Mohantv et alJl2002) . Since 
photometric variability can be caused by both magnetic activ- 
ity (through star spots) and dust clouds, it is not clear whether 
variability is more frequent or less frequent in M dwarfs than 
in L dwarfs. Measurements of v sin ; values show that later-type 
objects on average rotate much faste r and lead to expected ro- 
tation periods of about 1.5 - 13hr ( lBailer-Jonesii»200 4 ) for L 
dwarfs, whereas M dwarfs can have significantly longer pe- 
riods of up to two days. Thus it appears that the spin-down 
timescale is much longer for later-type objects. 

Stellar/substellar rotation periods can be measured directly 
from photometric time series if the objects show surface bright- 
ness inhomogeneities. This also permits a study of multiperi- 
odicity and/or nonperiodic variability. Candidates for the sur- 
face features are magnetically induced spots and, because of 
their lower effective temperature, dust clouds for late M and 
L dwarfs. By modelling the variability at different wavelengths, 
we can attempt to infer the physical properties of the atmo- 
sphere and surface features in individual cases. The few at- 
tempts to attribute a specific physical cause of the variability 
(via time-resolved spectroscopy) have been rather inconclusive 
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llBailer-Jonesll2002t IClarke et aljboojt iBailer- Jones & LamnJ 
l2003h . however, partly due to the high multifrequency photo- 
metric sensitivity (1 % or better) required to distinguish be- 
tween mechanisms. 

There have been several one- or two-band photometric sur- 
veys of ultra cool dwarfs (UCDs) of both young cluster objects 
and older field objects, mostly in the I-band but also in the R 
or near infrared bands. A comparison of these reveals variabil- 
ity (periodic or nonperiodic) to be present in about 40 % of 
objects in each of these two age groups (Bailer- Jones 2005). 
Interestingly, periodic variability is more common among the 
younger cluster sample, whereas cases of nonperiodic variabil- 
ity are more frequent in the older field sample. This could, in 
principle, be either an age or a spectral type effect. An age ef- 
fect might be related to a decline in activity, possibly from the 
dissipation of disks, although we have no explicit mechanism 
for this. A spectral type effect could be related to the onset of 
dust formation and, specifically, the dynamics of large scale 
dust clouds: Bailer- Jones & Mundt (2001, henceforth BJMOl) 
proposed a "masking hypothesis" in which the rapid formation 
and evolution of dust clouds masks the rotational modulation 
of the light curves of L dwar fs, thus accou nting for their non- 
periodic variability (also see iBailer-JonesI ((2004j)). We cannot 
distinguish between these two effects in the samples to date, 
because there is a broad age-spectral type (SpT) correlation in 
the object selection: most M dwarfs monitored were in young 
clusters (1-100 Myr) whereas the L dwarfs were field objects 
(probably more than several hundred Myr in age). Specifically, 
in our earlier work (BJMOl) - where we monitored 21 M6-L5 
dwarfs - we found variablity to be more common for objects 
later than M9. However, as stated in that paper, because of a 
coarse age-SpT correlation in the sample, we could not dis- 
tinguish between an age and a SpT dependence based only on 
those data. 

The main objective of the present work is to remove this 
correlation by extending the survey of BJMOl to (older) field 
M dwarfs. This represents a control sample against which the 
occurrence, periods, amplitudes and nature of variability in the 
field L dwarf sample of BJMOl can be compared, i.e. two pop- 
ulations with similar ages. Early and mid M dwarfs do not show 
photospheric dust formation (either theoretically or observa- 
tionally), so together with the L dwarf samples allows us to 
study the link between variability and dust. In addition to sig- 
nificantly increasing the time-resolved data on field M dwarfs 
in G, R and I, we are also able to constrain the physical cause 
of the variability in individual cases. 

The structure of the paper is as follows: The next section 
describes the target selection and data acquisition process. In 
Sectionl^we present the basic data reduction steps, photometry 
and a reliable error description model which together allow us 
to investigate low-amplitude variability. Methods for analysing 
general and periodic variability are discussed in Section 0] 
here we also describe our Monte Carlo approach to estimat- 
ing the reliability and uncertainty in the period detection pro- 
cess. We briefly mention the multichannel detection of a mas- 
sive flare in an M9 dwar f. This is reported in more detail in 
iRockenfeller et all (l2006l) . Section|5]deals with the overall re- 
sults plus details on some individual targets. Two of these ob- 
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Figure 1. The total efliciencies of the BUSCA G, R and I chan- 
nels, taking into account the dichroics, the CCD efficency and 
the Bessell-I filter used in the I-band. Also plotte d is the spec- 
trum of the M4.5 dwarf DENIS P-Jl 158-1201 of iMartm et all 
( 19991 

jects are studied for the cause of the observed variability and a 
comparison of the frequency of variability for L and M dwarfs 
is performed. We finish with the conclusions in Section |6l 
For a more detailed description and discussion o f the methods 
and for comments on every individual target see iRockenfelleJ 
( 20051 

2. Data Acquisition 

2.1. Target selection 

The target list was assembled by selecting bright M dwarfs 
{mi < 17.0 mag) which are visible from the Calar Alto ob- 
servatory, Spain, in June for several hours at an airmass of less 
than 2.0. Many of the objects are bright enough such that the 
G-band data have a relatively high signal-to-noise ratio (SNR). 
Strong activity as measured by was avoided and they were 
chosen to cover most of the spectral M-type sequence (M2 to 
M9). A total of 19 field M-dwarfs were observed over the two 
observing seasons. Target details are given in Tabled] 

2.2. Observations 

To perform simultaneous multiband photometry, we used the 
BUSCA four channel CCD camera at the 2.2 m telescope at 
Calar Alto Observatory, Spain. This instrument uses dichroics 
to split the light beam into four wavelength bands, namely UV, 
G, R and I. The UV, G and R passbands are defined by the CCD 
response and dichroic transmission function; for the I band we 
additionally used a Bessel I filter (see Fig. [Q. Each channel 
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Table 1. Target properties: full name, spectral type, I magnitude, number of images in total, year of observation, nights of 
observations, duration of the observations on the individual nights, additional information and the references to these. Objects 
are sorted by increasing right-ascension. 



target 


SpT 


mi 


Nobs 


year 


nights 


duration [hr] 


further information 


reference 


2MASSW J1311391+803222 


M8 


15.9 


33 


2002 


5,8 


5.9, 6.7 


H„ = 3.0 


GOO 


2MASSW J1336504+475131 


M7 


15 


66 


2003 


9,10 


4.6, 4.8 


V sin / = 30, = 3.2,5.0 


GOO, R02 


2MASSW J1344582+771551 


M7 


15.2 


41 


2002 


1,4 


5.2, 6.2 


H„ = 2.7 


GOO 


LHS370 


M5 


12.3 


20 


2002 


8 


6.1 






LHS2930 


M6.5 


13.3 


58 


2003 


7,8 


7.1, 7.0 


R/Rsun ~ 0.33, Teff = 2687 ± 65 K 


H04, D02 


CTI153948. 1+280322 


M2 


14.9 


50 


2003 


1,3 


4.9, 6.9 


Binary Star; v sin / = 7.0, = 0.7 


M03 


2MASSW J 1546054+374946 


M7.5 


15.2 


51 


2003 


4,5,6 


6.4, 2.6, 2.4 


H„ = 10.9 


GOO, C03 


LHS3189 


M5.5 


14.9 


49 


2002 


1,2 


6.5, 6.7 




B91 


2MASS J 16272794+ 8 105075 


M9 


16.4 


52 


2003 


4,5,6 


6.4, 2.9, 2.4 




GOO 


CTI162920.5+280239 


M4 


15.6 


49 


2003 


1,3 


5.2, 6.4 


H„ = 10.9 




2MASSW J1707183+643933 


M9 


15.9 


34 


2002 


4,5 


6.4, 5.6 


T,tf ~ 2350K, //„ = 9.8, 
logg = 5.16, M/Mi,,,, = 0.062 


GOO, G03 


CTI170958.5+275905 


M5.5 


13.9 


54 


2003 


4,5,6 


6.7, 2.8, 2.3 




K94 


2MASSW J1714523+301941 


M6.5 


14.9 


51 


2003 


1,3 


6.1,5.4 


Teff a 2775 K, V • sin ; < 4, 
H„ = 3.2, 5.4 


GOO, R02 


LHS3307 


M5 


15 


57 


2003 


9,10 


6.8, 4.1 




B91 


2MASS J17501291+4424043 


M7.5 


15.6 


64 


2003 


7,8 


6.8, 6.9 


H„ = 2.7 


GOO 


LHS3339 


M6 


14.0 


50 


2002 


1,2 


6.0, 6.4 


Teff = 2957 ± 70 K, R/Rs,„ « 0.33 


LOO, D02 


2MASSW J 1757 154+704201 


M7.5 


14.2 


38 


2002 


4,5 


6.3, 7.0 


= 3.0 


GOO 


CTI180120.1+280410 


M2.5 


14.8 


57 


2003 


7,8 


6.2, 6.9 






LHS3376 


M4 


10.7 


30 


2002 


8 


6.0 


Teff ~ 3100 K, V • sin J = 14.6 ± 1.0 


D98 



B91: Bessell Q99r); K94: Kirkoatrick et al. (1994); D98: Delfosse et al. (1998) ; GOO: Gizis et al. (2000) 
LOO: Leggett et al. (2000 ): R02: Reid et al. (2002); D02: Dah n et al. (2002) 

M03: Mohantv & Basri (2003) : C03: Cruz et al. (2003) ; G03: Gorlova et al. (2003) : H04: Henry et al. (2004) 



is equipped with a CCD485 Lockheed Martin 4Kx4K CCD. 
The CCDs are thick, except for the UV channel which uses a 
backside-illuminated thinned chip. For all CCDs 2x2 binning 
was used, reducing read-out time and leading to a pixel scale 
of 0.352". The field-of-view of 12' ■ 12' allowed us to select 
a large number of reference stars for differential photometry. 
For use throughout this paper, we define an adjusted Julian day 
(AJD) by AJD = JD - 2 450 000. Observations were conducted 
in two runs: the first with 8 usable nights in June 2002 (AJD 
2424.4-2431.7) and the second with 9 usable nights in June 
2003 (AJD 2794.4-2803.6). On each night, two or three tar- 
gets were observed alternately and for each target, data were 
collected on two or three nights. Integration times of 45 to 300 s 
were used depending on target brightness and weather condi- 
tions, to achieve high signal-to-noise data and still retain sensi- 
tivity to periods in the sub-hour regime. 

3. Data reduction 

3.1. Basic reduction steps 

All science images were reduced using the IRAF' package. 
UV-band data could not be used due to the faintness of the 
M dwarfs in this spectral range, other than for finding flares 
in the targets' time series (see Section R3t . We reduced the 
other three channels with the same reduction steps, with pa- 
rameters set to achieve homogeneous, comparable photometry. 

' The Image Reduction and Analysis Facility, provided by the 
National Optical Astronomy Observatories (NOAO). 



The basic reduction steps include an overscan subtraction and 
flat fielding process. Because non-thinned CCDs were used, no 
fringes can be seen in any of the images and hence fringe cor- 
rection was not performed. 

The zero integration time frames, which were obtained to 
correct for two-dimensional bias patterns in the science images, 
are useless: They all show vertical bars caused by the read-out 
electronics which do not appear in the science frames, a known 
problem that occurs frequently when using BUSCA. To partly 
compensate for this, a first-order fit to the overscan area was 
subtracted from each science and flat field frame. This corrects 
at least for one-dimensional patterns in the y-direction of the 
bias level. Because residual two-dimensional patterns are only 
weak and do not show a strong gradient on the scale of the ap- 
plied sky apertures, their contribution to the total photometric 
eiTor is well below 0.5 %. We used an clipped average of ten 
dome flat fields to correct for changes in the pixels' quantum 
efficiencies. Individual flat field frames were applied separately 
for each channel and each night. 

3.2. Pliotometry 

To achieve the desired accuracy of better than 0.01 mag, dif- 
ferential photometry was performed to reduce the effects of 
temporal variations in Earth's atmosphere. Aperture photom- 
etry was done on the target and tens of reference stars; these 
were chosen according to the criteria used in BJMOl. Refer 
also to that paper for a description of differential photome- 
try. Although various aperture radii were tested, the results we 
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show are with a radius of 6 pixels (except for LHS370, see 
Section which we found to maximize the signal-to-noise 
ratio (SNR). The large field-of-view guarantees that in most 
cases more than 15 (and even up to 50) non-variable reference 
stars could be found within the target field, ensuring a stable 
reference flux. For the procedure of creating the reference star 
list, see the next section. We excluded the presence of sign ifi- 
cant se cond-order-extinction (SOE; see lSailer-Jones & Lamml 
( I2OO3I) for a discussion) in our data by: plotting relative mag- 
nitudes versus airmass; comparing the relative light curves 
formed with only blue and only red reference stars; plotting 
the variability measure of all reference stars of a field (y^ val- 
ues, se e Section lTTl versus the R - I colour. See .RockenfelleJ 
(I2OO5I) for more details. 



3.3. Error sources and estimation 



An accurate error estimate is important for assessing the 
presen ce of variability. As discussed in Bailer-Jpnes & Mund^ 
lll999l) . the total photometric error is composed of formal 
(Poisson noise) and informal error contributions. As imperfect 
flat fielding is thought to dominate the informal error, we first 
tried to model the total photometric error (o-tot) by adding a 
constant term (of 0.5 %) in quadrature to the theoretical (for- 
mal) errors provided by IRAF (itiraf)- However, this lead to 
an undesirable magnitude dependence of the variability mesi- 
smex^- Instead (and to overcome this), we finally used the fol- 
lowing error model: because all reference stars of a target field 
are non-variable, the scatter (standard deviation cTims) in their 
light curves is a measure of the total photometric error at the 
corresponding magnitude. We plotted 0-^^^ versus (Tiraf and 
fitted a first order polynomial to this graph. Applying this fit, 
i.e. (Ttot - a + b ■ (Tiraf where a and b are the fit's free parame- 
ters, we arrive at a reliable error description. Typical values of a 
and b are -0.01 to 0.01 and 0.5 to 1.5, respectively. The fitting 
procedure was done for each field and channel separately. The 
reference star selection is an iterative process, which means that 
first the;^'^ measure (see Section im is determined for all can- 
didates. Then the most variable one is excluded from the set and 
the process is repeated until no variable star remains. Our er- 
ror model partly accounts for varying data quality on different 
nights, e.g. due to bad seeing conditions, since this manifests 
itself as an increased standard deviation in the reference stars. 



We further found that other informal errors contribute far 
less than 0.5 %. This includes imperfect overscan subtraction 
and two rather unusual sources: weak charge trailing on the 
CCD for brighter stars (and which appears only in the G- 
channel) and what we call glitches. These are quite similar in 
size and intensity to cosmics but are caused by the read-out 
electronics. They are quite numerous in the I-channel (some 
hundred per image) but only rare in R and G. 



4. Time series analysis 

4.1. test 

To test whether or not deviations in a target's relative light 
curve are consistent with the photometric errors (the null hy- 
pothesis) we evaluated the measure 




where mrei(0 is the relative magnitude in the /th (of A^) frame 
and 5mrei(0 the error therein. The larger the;if^ value, the larger 
the probability that the null hypothesis is wrong and that the 
object is variable. We claim an object to be variable if the prob- 
ability for the null hypothesis p is smaller than 0.01. We used 
this test first to assemble the set of reference stars and then to 
test for possible variability within the targets. 

We often encountered the case that significant general or 
periodic variability is only present in one or two channels. This 
can either be due to different sensitivity limits or to different 
variability amplitudes in the channels (or to a combination of 
both). As we will discuss in chapter|5l surface features can in- 
deed lead to such a behaviour. 

4.2. Periodic variability 

Periodic variability can be caused by co-rotating surface fea- 
tures which are stable on time scales of an object's rotation 
period. To check for periodic behaviour, the following scheme 
was applied to each target: the CLEAN periodogram was cal- 
culated and searched for significant peaks; in case there is such 
a peak, the target's light curve was phased to the corresponding 
period and checked if it confirms the period in question. The 
Lomb-Scargle periodogram was also evaluated for all targets 
but the influence of the spectral window function is very strong 
for periods longer than approximately 10 hours. Because of 
this the Scargle power spectra are not used for the final results. 
Nevertheless, the shorter periods of 2M1707+64 and LHS370 
are confirmed by this method. 

The CLEAN algorithm tries to remove the influence of the 
discrete and finite sampling of observational data on power 
spectra . For more information on this, see BJMOl and Ro bert si 
( Il987t) . To judge whether a peak in the power spectrum is sta- 
tistically significant, we performed Monte-Carlo simulations 
to determine the false alarm probability (FAP) power levels 
for each time sampling. The FAP denotes the probability that 
a peak in the power spectrum is caused by nois e. We chose 
to use the same methods as already discussed bv iLamm et alJ 
Artificial light curves of non-variable stars were created 
by simulating pure Gaussian noise as well as by shuffling the 
actual magnitudes of the target's light curve with respect to the 
real epochs. For both methods, the highest peak in the corre- 
sponding CLEAN periodogram was determined and the power 
level that was exceeded by 100 out of the 10"^ simulations was 
defined to be the 1 % FAP power level (similarly one obtains 
the 10 % and 0.1 % levels). To be conservative we compare the 
corresponding two values obtained for both methods and use 
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Figure 2. Period uncertainty in hours plotted versus the loga- 
rithm of the period in hours; obtained for an amplitude ratio of 
1.5 for the 2M1707+64 I time sampling with the method de- 
scribed in section lT^ (solid line and solid circles). The results 
of equation[2are also included as open circles on a dashed line. 



the larger one. We claim periodic variability if peaks above the 
1 % limit are present. 



tions is that they assume a sinusoidal signal, whereas periodic 
variability does not necessarily match this. 

A similar procedure to the above has been used by 
IScholz & EislofFel ( 2004a b) who find a similar behaviour of 
the period uncertainty. Also shown in Fig. |2] (open circles) 
are the results of a theoretical period uncertainty estimate (for 
small eiTors): 



Here Av is the width of the main peak of the window function 
W(v), that can be approximated for not too uneven data sam- 
pling by Av « J. Here, T is the total time span covered by the 
observations. P and AP are the period in question and its uncer- 
tainty, respectively. For long periods (larger errors) the approx- 
imation made to derive the above equation may no longer be 
valid. This would explain the huge difference relative to the re- 
sults obtained with the simulations at longer periods. For more 
information on equation^ see Roberts ( 1987). 

Since the predicted uncertainties may differ by a factor of 
up to 3 between these two methods and since both methods 
do have drawbacks, it is difficult to decide which one to use. 
Thus the period uncertainties we state are always the results of 
the simulations but we also give the values obtained by equa- 
tion^in parenthesis. Apart from just finding errors on the peri- 
ods found with this survey, the period uncertainty is important 
when periods are reported in more than one channel of a tar- 
get. If those periods are close to each other, it is likely that they 
actually correspond to a single periodicity, e.g. the rotation pe- 
riod, if their values do not differ by more than two or three 
times their uncertainty. 



4.2.1. Period uncertainty 

To investigate the uncertainty in the detected periods, we sim- 
ulated a sinusoidal signal at random phase with respect to 
the time sampling and added Gaussian noise. We define the 
Amplitude Ratio (AR) to be the ratio between the root-mea- 
square (rms) amplitude of the sinusoidal and the one of the 
noise. The absolute value of the difference between the period 
of the input signal and that of the highest peak in the CLEAN 
power spectrum was averaged over 10'* simulations at each in- 
put period. This result is a measure for the expected period 
uncertainty at the period of the sinusoid. An example for one 
time sampling is shown in Fig.|2l The input frequency is varied 
over the whole range of frequencies (actually we divided that 
range into 100 bins and performed lO'* simulations per bin). 
The shape of the curve is only weakly dependent on the AR 
for values above 1 .5 and period uncertainties only slightly in- 
crease for ARs down to 1.0. Table0]shows that, except for one 
case, all claimed periods occur at ARs of more than one. The 
simulations show that the period accuracy is very good for pe- 
riods up to about 10 hours and becomes much worse between 
15 and 30 hours. This can be attributed to the gap in the data 
between the two consecutive observation nights. Although we 
expect the uncertainty to become lower for yet longer periods 
(since we observed the target on two nights), it is surprising 
how low it actually gets. A possible problem of these simula- 



4.2.2. Detection Fraction 

To study the sensitivity of our period detection procedure, we 
introduce the detection fraction (DF). This quantity is equal to 
the fraction of detected periodic signals as a function of period; 
it is calculated using simulations of a sinusoidal and Gaussian 
noise (with the same definition of the AR as in the last section). 
The whole range of frequencies was divided in 100 bins and the 
simulations performed separately for each bin. The fraction of 
simulations with a fixed input frequency which lead to a peak 
in the power spectrum above the 1 % FAP power level is equal 
to the DF. We applied two different DFs, one where peaks any- 
where in the periodogram were considered (DFaii) and another 
where a significant peak was counted only if it was found inside 
the input bin (DFbin). Naturally, the detection fraction depends 
strongly on the AR. See Fig.|3lfor the results on the same target 
as in Fig. 121 The resemblance to the plot of the period uncer- 
tainty is striking and naturally because both quantities depend 
on the same data sampling. Thus the same reasons can be given 
for the minimum in the DF at about f = 10 - 30 hr and for the 
maximum in the period uncertainty at the same period range. 
The difference between the solid and the dashed line in Fig.|3lis 
the fraction of period detections at the wrong frequencies, i.e. 
outside the input bin. In agreement with the period uncertainty 
simulations, this difference is quite small at periods shorter than 
10 hours and thus if we detect a periodicity within this range 
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0.5 1 
logio(period [hr]) 

Figures. The plot shows various detection fractions for the 
2M1707+64 time sampling over the logarithm of the period in 
hours. Solid lines represent detection fractions which require 
a peak to be found in a narrow range around the input period 
(i.e. inside the input bin), (DFbin), and dashed lines the ones 
that allow a peak to be at any period (DFaii). For periods of 
up to logio(period) =1.0 three amplitude ratios are investi- 
gated: 3, 2, 1.5 from top to bottom within the plot; for each 
ratio one solid and one dashed line. Because the general be- 
haviour does not change, and to avoid crowding, for periods 
longer than ten hours only the curves corresponding to an am- 
plitude ratio of 2 are shown. Here the solid and dashed lines 
merge at logio(period) = 1.47. 



it is very likely detected at the right period. If we again com- 
pare the ARs given in Table0]with the Detection Fraction plots 
of the individual targets, we can summarise the results as fol- 
lows: for most time samplings, an AR of 1.5 leads to a DF 
between 50 % and 80 %; for ARs of 1.0, the DF Ues between 
15 % and 30 % (in the case of 2M1707H-64, see Fig.Q the DFs 
are somewhat lower at each AR). Hence it is not unlikely that 
we missed periods for targets with ARs of less than 1.5, even 
in the range of up to ten hours. Furthermore, the vsin; dis- 
tribut ion of M dwarfs (Mohanty & Basri 2003t iDelfosse et alJ 
Il998l) suggests that there are some objects, particularly at early 
M type, with expected periods of much more than 10 hours (up 
to a few days). For these, the sensitivity of our detection proce- 
dure is questionable. Hence, to lower the probability of missing 
periods significantly, one would have to perform higher SNR 
observations and/or get a longer time base. 



4.3. Flares 

We will very briefly address the topic of flares in the observed 
field M dwarfs. Fo r detailed information on t he flares detected 
in this data set, see ' RockenfeUer et alJ tOOdi . 

Although well-known, the exact processes taking place dur- 
ing a flare event are not well-understood. Probably, magnetic 
energy is transfered to thermal energy and thus leads to a 
brightening of the effected area on or near the stellar surface. 
Flares on the Sun are often associated with eruptive promini- 
nences, sometimes ejecting solar material and charged particles 
into the solar system. 

In principle, multiband monitoring in optical bands is ideal 
to detect flares because the amplitude of the brightness vari- 
ations increases tremendously from the I- to the UV-band. 
A flare is characterised by a fast rising signal followed by a 
slower decreasing one (of exponential shape). The duration of 
the event is positively correlated with its amplitude and ranges 
from a few minutes (or even shorter) to a few hours. This cre- 
ates the problem that low-amplitude flares are too short to be 
seen in detail in this data set since the minimum time span be- 
tween two data points is about 5 minutes. However, one huge 
flare was detected in 2M 1707 +64, with an UV-amplitude of 
more than 6 magnitudes and a recorded duration of about 1 
hour For this event we captured the brightness evolution over 
five data points. Three other events that are probably flares were 
found in 2M1714+30, 2M1546+37 and 2M1344+77 at lower 
amplitudes (between 1.5 and 2.7 mag in UV). The total obser- 
vation time of this survey is 218.0 hours, yielding a flare rate of 
9.2- lO^'^hr"', 1.83- lO"-' hr"' if we only count the two strongest 
events or all four, respectively. 

5. Results 

5.1. General results 

The results of the test and the period search for all 19 tar- 
gets are listed in Tabled It shows the spectral type, the number 
of finally used reference stars (A^,s), the value of the light 
curve and the one corresponding to p = 0.01 (Xq), variabil- 
ity flags, the probability of the null hypothesis ptot, the vari- 
ability amplitude, a possible period, the Amplitude Ratio and 
the finally claimed variability flag. Here, the variability ampli- 
tude is simply the root-mean-square (rms) of the relative light 
curve in case of variable objects and an estimation of an up- 
per limit above which variability would have been detected for 
non- variable objects (using the method of BJMOl, Section5.1). 
The variability flags in column 6 state whether or not a target is 
variable according to the x^ test on the individual observation 
nights (and in total in parentheses). In contrast to this, the flag 
in the last column indicates whether we finally claim the target 
to be variable, after consulting all means of investigation. We 
claim a target to be generally variable if the probability of the 
null hypothesis of the x^ test is less than 1 %. This choice is 
quite arbitrary and using other limits would turn some detec- 
tions into non-detections or vice versa. Likewise, periodic vari- 
ability is claimed to be significant if the corresponding peak 
in the periodogram has a FAP of less than 1 %. In all cases 
where we found variabihty or remarkable features in a target's 
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Figure 4. Relative light curves of LHS370 at an aperture radius 
of nine pixels. The I, R and G band time series are shown with 
a solid, dashed and dotted line, respectively. Typical error bars 
are plotted for each channel. 



light curves, we also visually checked the reference stars' light 
curves, variability measures and power spectra to rule out ex- 
terior influences. 

5.2. Discussion of individual targets 
LHS370 (M5) 

Due to bad seeing of up to 3" on night 8 in 2002, the star aper- 
ture radius had to be increased for LHS370 from 6 to 9 pixels to 
maximise the SNR. This object shows signs of variability in all 
three channels. In I it is a strong detection (p < 10"*) accord- 
ing to the test but a non-detection in R and G. Nevertheless, 
the light curves (Fig. |3 of R and G show variations at sim- 
ilar amplitudes to that in I, but with larger eiTors. Although 
our data only cover about one cycle, the shape of the light 
curves suggests a periodic behaviour except for the I-band. 
This is confirmed since the power spectra do contain signif- 
icant (above the 0.1 % FAP power level) peaks in R and G 
at periods of 5.9 + 2.0 hr (theoretical uncertainty 2.9 hr) and 
6.5 + 2.0 hr (±3.5 hr), respectively, but none in I. See Fig.|5]for 
the CLEAN periodogram of the G-band data and its light curve 
phased to the just mentioned period in Fig.|6l Since the phased 
light curves in both R and G look reasonable, there is good 
evidence for periodic variability. Studying the results of the pe- 
riod uncertainty simulations (see above), we find that they are 
consistent with a single periodicity. Since the Amplitude Ratio 
(AR) of the I-band data is 2.4 we can confidently infer that 
there is no similar period in the I-band compared to those in 
R and G. Because of the higher AR in R (compared to G), we 
claim the R-band period at 5.9 ± 2.0 hr to be the object's rota- 
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Figure 5. CLEAN periodogram of the G-band data of LHS370. 
The upper panel shows the spectral window function and the 
lower panel the dirty power spectrum as a dotted line as well 
as the cleaned power spectrum as a thick solid line. The 0.1 %, 
1 % and 10 % FAP power levels are plotted as dotted horizontal 
lines (from top to bottom). 



tion period. This i s consistent with v sin ; rn easu rements of var- 
ious M dwarfs bv'Moha ntv & Basril (l2003l) and'Pelfosse e t alJ 
fl998). To investigate possible sources for this kind of vari- 
ability, we compared the observations with a grid of synthetic 
atmosphere spectra. For details on this procedure see the sec- 
tion on 2M1707+64 below. From such an analysis we could 
not match the observed variability amplitudes simultaneously 
in the three channels with any reasonable atmospheric model 
spectra featuring either magnetically induced star spots or dust 
clouds. Alternative models for dust clouds or spots may be re- 
quired. We define the coverage factor e as the fraction of the 
visible hemisphere that is covered by surface features. By sim- 
ply comparing the predicted amplitudes with the observed, we 
can however place an upper limit of e = 0.05 on the coverage 
factor of spots on the surface of LHS370. Otherwise we would 
have been able to clearly detect their signature in the I-band. 
Higher SNR data would yield a higher photometric precision 
and thus allow us to better constrain the source of variability in 
LHS370. 



2M1 707+64 (M9) 

This target shows various kinds of variability. First, it is a 
strong detection in terms of the test in I and R on both ob- 
servation nights {p < 10"* and p = 0.001, respectively) and 
only on night two in G. Besides this general variability indica- 
tion, periods with FAPs better than 0.1 % at 3.65 hr, 3.7 hr and 
3.3 hr are present in I, R and G respectively. The uncertainty 
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Figure 6. G band light curve of LHS370 phased to a period of 
6.5 hours. The typical error is also shown. Note that two cycles 
of the phased light curve are plotted. 

in those periods is estimated to be about 0. 1 hr (theoretical es- 
timate 0.6 hr) and thus all three values correspond to the same 
periodicity. As with LHS370, this period is consistent with the 
V sin / values of late M type objects, so we claim this to be the 
rotation period. Figure shows that these periods can be seen 
visually as nearly sinusoidal modulations in the light curves. 
Figure|8lshows the CLEAN power spectrum in I and Fig.Othe 
I band light curve phased to the mentioned period. A huge flare 
was detected at the end of the second night which confirms that 
magnetic activity is significant in this late type object. For more 
inf ormation on this and other flares found within this data set, 
see iRockenfeller et alJ ( l2006l) . 

Since the variations in the light curves of the different 
channels are strongly correlated, it is obvious to assume that 
some surface feature co-rotating with the target is the rea- 
son for the measured period. To put stronger constraints on 
the source of variability in 2M1707-H64, we u sed the syn- 
thetic atmosp heric models of lAUard et alJ(l200l V As shown in 
lBailer-.Tonesl (t2002) . it is straight forward to derive the differ- 
ence spectrum caused by co-rotating or forming / dissolving 
star spots or dust clouds. By integrating the predicted spectral 
variations over our photometric bands, we arrive at variability 
amplitudes for a specific model (refer to the just mentioned pa- 
per for a description of the two cases used, dusty and cond, as 
well as example spectra plots). Since, for a specific model, the 
predicted amplitudes depend linearly on the coverage factor e 
of the surface features (for values of up to e = 0.3), the ratio of 
two amplitudes in different wavebands will be independent of 
e. This makes these ratios suitable for assessing whether or not 
an individual model fits our data. Based on this, we built a grid 
of results for various trial atmospheric models. The best-fitting 
model (dusty atmosphere with T^g = 2300 K, logg = 5.0 and 



Figure?. The relative light curves of 2M 1707-1-64 in I (sohd 
circles on a solid line), R (open squares on a dashed line) and 
G (solid triangles on a dotted line). Each observation night is 
shown in an individual panel with typical error bars for each 
channel. 



t i 


— ' ' 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _| 











0,5 1 1,5 2 2,5 

frequency [l/hr] 



Figures. CLEAN periodogram of the I band light curve of 
2M1707+64. The upper panels shows the spectral window 
function and the lower panel the dirty power spectrum as a dot- 
ted line as well as the cleaned power spectrum as a thick solid 
line. The 0.1 %, 1 % and 10% FAP power levels are plotted 
as dotted horizontal lines (from top to bottom). Cf. Fig.|5] The 
above periodogram has higher resolution due to the longer time 
base. 
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Table 2. Comparison of observed and simulated variability amplitudes. The upper panel shows the observational standard (std), 
i.e. RMS, and peak-to-peak (ptp) amplitudes and the ratios formed hereof. Ratio errors are estimated from the errors in the 
relative magnitudes. Theoretical amplitudes and ratios are presented in the lower part. The four cases are cool spot or dusty cloud 
on a cond (clear) atmospheres, and cool spot or a hole (cond "cloud") on a dusty atmosphere. The underlying atmosphere has 
an effective temperature of 2300 K and surface gravity of \ogg = 5.0; the cool spot is 100 K cooler and all feature types (spot or 
cloud) have a projected surface coverage factor of 0.1. These are the parameters of the best-fitting model. Values in parentheses 
indicate the amplitude ratios as explained in the last column. Note that the ampUtude ratios are independent of the coverage factor 
for small (< 0.3) coverage factors. 







observational quantities 






band / colour 


std ampl. [mag] 


ptp ampl. [mag] 


ratio 


std ratio 


ptp ratio ratio error 


G 


0.022 


0.055 


G:R 


1.57 


1.49 ±0.25 


R 


0.014 


0.037 


G : I 


1.83 


1.67 ±0.25 


I 


0.012 


0.033 


R : I 


1.17 


1.12 ±0.30 






theoretical quantities 








cond atmosphere 


dusty a 


tmosphere 


legend 


band / colour 


cool spot 


dusty cloud 


cool spot 


cond cloud 


ampl. (ratio) 


G 


0.039 (1.56) 


0.064(1.42) 


0.072 (1.48) 


-0.135 (1.85) 


G(G:R) 


R 


0.025 (1.56) 


0.045 (3.05) 


0.049 (1.70) 


-0.073 (5.40) 


R (G : /) 


I 


0.025 (1.00) 


0.021 (2.14) 


0.042 (1.15) 


-0.025 (2.92) 


l(R:l) 



Table 3. Supplementary data (to Tabled) that shows the amplitude ratios as calculated with model parameters given in the first 
two columns (the surface gravity is fixed at \ogg = 5.0 since it does not have a significant influence on the results; feature sizes 
are independent of the coverage factor up to e 0.3). cond spot stands for a cool cond spot on a cond background atmosphere 
and dusty cloud similarly denotes a dusty cloud on a cond atmosphere. 



grid of theoretical quantities - amplitude ratios 
model parameters cond atmosphere dusty atmosphere 





(Ar),p„,[/f] 




G:R 


G : / 


R : 1 




G:R 


G : / 


R : I 


2300 


100 


cond spot 


1.56 


1.56 


1.00 


dusty spot 


1.48 


1.70 


1.15 






dusty cloud 


1.42 


3.05 


2.14 


cond cloud 


1.85 


5.40 


2.92 


2300 


200 


cond spot 


1.62 


1.38 


0.85 


dusty spot 


1.24 


1.33 


1.07 


2300 


300 


cond spot 


1.60 


1.24 


0.77 


dusty spot 


1.14 


1.16 


1.02 


2500 


100 


cond spot 


1.22 


1.92 


1.57 


dusty spot 


1.27 


1.96 


1.54 






dusty cloud 


1.44 


6.82 


4.73 


cond cloud 


1.59 


8.70 


5.46 


2500 


200 


cond spot 


1.22 


1.70 


1.39 


dusty spot 


1.21 


1.66 


1.38 


2500 


300 


cond spot 


1.22 


1.49 


1.22 


dusty spot 


1.15 


1.38 


1.20 


2100 


100 


cond spot 


2.25 


1.16 


0.51 


dusty spot 


1.36 


1.23 


0.91 






dusty cloud 


1.21 


1.60 


1.33 


cond cloud 


2.57 


5.50 


2.14 



a lOOK cooler spot) along with the observational amplitudes 
are shown in Table |2j whereas Table |3l shows some other re- 
sults from the computed model grid. The theoretical predic- 
tions were made for cooler spots on either type of background 
atmosphere (cond and dusty) as well as for clouds of the op- 
posite type and the same temperature as the background atmo- 
sphere, i.e. a cond cloud on a dusty atmosphere (clear hole in 
a dusty sky) and the other way around (dusty cloud on a clear 
sky). This setup excludes fast convection as a source for the 
clouds since they would then be hotter than the surrounding re- 
gions. Because of their small size (compared to the noise in in- 
dividual measurements), the determination of the observational 
amplitudes is non-trivial. We decided to form both the peak- 
to-peak amplitudes (ptp) after an exclusion of obvious outliers 
and the root-mean-square of the light curve, here called stan- 
dard amplitude (std). In the case of 2M1707+64, the amplitude 



ratios formed with std and ptp amplitudes are similar because 
of the nearly sinusoidal shape of the light curves. For a per- 
fectly sinusoidal signal the std amphtude is y/n/l times the ptp 
amplitude. 

Studying both tables for T^ff = 2300 K and 7,^ = 2500 K, 
we can clearly exclude clouds as the source of variability be- 
cause the predicted signatures are too far from the measured 
ones. The case of a dusty cloud on a T^ff = 2100 K cond atmo- 
sphere also is an acceptable fit to the data. But first, the effec- 
tive temperatur e is somew hat lower than the one derived by the 
spectral type ( Gorlova et al..(.2003) . although the uncertainty in 
Teft is about 200 K) and second, and more convincingly, cond 
atmospheres are expected to occur not earlier than T-type ob- 
jects. 

Besides the best-fitting one, models with spots on (particu- 
larly dusty) atmospheres of 2200 < T^ff < 2500 fit the observed 
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Figure 9. Relative light curve of 2M1707+64 in I phased to a 
period of 3.65 hr. The plot shows the combined data of both 
nights. The typical error is also plotted. Note that two cycles of 
the phased light curve are plotted. 
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Figure 10. Relative Hght curves of CTI1709+27 in the 1-band 
(solid line), R-band (dashed line) and G-band (dotted line). 
Each observation night is shown in an individual panel with 
typical error bars for each channel. 



amplitude ratios within the estimated errors. Hence we con- 
clude that co-rotating magnetically induced spots, most likely 
on a dusty atmosphere of 2200 - 2500 K, cause the observed 
variability in 2M1707-H64. We infer from the observed ampli- 
tudes that the coverage factor of these spots cannot be larger 
than about e = 0.075 for a dusty atmosphere (and about 0.15 
for a cond one). Note that a coverage factor of e.g. e = 0.1 does 
not necessarily refer to a single spot covering ten percent of the 
visible hemisphere. One can imagine other spot configurations 
leading to the same spectral signature as that of a single spot. 
For example, a symmetrical distribution of many small spots 
(with the same e) will not result in any detectable photomet- 
ric signal. Thus the value of e = 0.075 found for 2M1707+64 
could also mean that a larger fraction of the star's hemisphere 
is covered by symmetrically distributed spots and that there is 
a non-symmetric spot coverage of about 7.5 %. 

cm 709+27 (M5.5) 

We present the relative light curves of this target in Fig. 
as an example for a typical object showing general variability 
above a probability of 99 % (here even at p < 10"*, see Table|^ 
according to the;i^^ test. No remarkable features or periodicities 
are present. 

The choice of the probability at which to claim a target to 
be variable in terms of the test is quite arbitrary. The use of 
p - 0.05 instead of p - 0.01 would have led to 7 more targets 
to be claimed variable (at least in one channel each)! Thus one 
should first use a fairly low value for p and second, it would 
be ideal to re-monitor those targets after some time to test for 



persistent variabiUty. Of course, surface features could appear 
only temporarily and thus lead to a non-detection during the 
re-monitoring, even though they actually were present on the 
first observation run. Finally, this problem can be overcome by 
higher SNR data that permit even lower amplitude detections. 

Further variable objects 

Hints of periodic variability were also found in 2M 1344+77 
in G (12.5 + 3.0 hr), in LHS2930 in 1 (13.2 + 1.9 hr) and 
in 2M1714+30 in R (6.9 + 0.15 hr) by the evaluation of the 
CLEAN power spectra, although the the visual appearance of 
the phased light curves do not support the presence of such pe- 
riods. On the other hand, simulations of periodic modulations 
show that we should not always expect this, especially with 
multiple spots: see section 6.2 of BJMOl. Hence Table 0]lists 
these periods with a question mark as these periods are some- 
what tentative. Of these three targets only 2M1714+30 (in R) 
is a detection according to the test. 

General variability according to the test was detected 
in other targets as well, but they do not show any remarkable 
features. These targets are: CT11539+28 (1), 2M1546+37 (R), 
CT11629+28 (R and G), LHS3307 (R), LHS3339 (I and G), 
LHS3376 (1). VariabiUty in the G-band data of CTI1629+28 
was rejected because the two targets which were also observed 
in the same nights show similar trends in their corresponding 
G-band light curves. Because of bad data quality in the case of 
LHS3376, the variations in 1 also cannot be considered to be 
significant. 
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Significant correlations between different channels could 
only be found for 2M1707+64 and LHS370, where they are 
obvious. Correlations are another indicator of the reliability of 
a variability. Of course, the underlying physical process could 
produce polychromatic variations such that the signal is too 
small to be detected in one or more channels. 

In total, we find non-periodic variability in 5 out of 19 
M dwarfs (or 4, if we consider the period of 2M1714+30 to 
be real). Periodic variability was convincingly found in 2 tar- 
gets (or up to 5, if we choose to be less conservative). Although 
these results are surely affected by small number statistics, they 
suggest that about 2/7 ^ 30% of variable field M dwarfs show 
periodic variability. 

The fact that we did not detect any variability for 12 out of 
19 targets means that either no detectable surface features are 
present on these objects or that their rotation periods are too 
long to be detected. The latter may be particularly relevant for 
the four earliest type objects in our sample with spectral types 
M2-M4, since rotation periods of up to a few days are possi- 
ble. The large fraction of non-periodic variable targets (5 out 
of 7 in total) could be due to rapidly evolving surface features 
which would veil the rotation period (the "masking hypothesis" 
ofBJMOl). 

It is interesting to note that in the case of LHS370, no peri- 
odic variability was found in the 1-band even though we assign 
a rotation period to this target on account of the R- and G- 
band data. Hence we would have counted LHS370 to be gen- 
erally variable but non-periodic if it were only observed in the 
I-band. Since almost all recent monitoring programs on UCDs 
have been conducted only in the 1-band, it is possible that the 
large number of non-periodic variable L dwarfs found would 
show periods at other wavelengths. Multichannel observations 
of L dwarfs are required to investigate this. 

Non-variable objects 

All other targets show no signs of variability or exceptional 
behaviour and are hence not mentioned individually. At the 
present sensitivity Umit, we consider them non-variable. See 
Tablel^for upper limits on variability amplitudes and for ARs. 
The latter allows us to judge whether we could have missed 
existing periods because of this, see Section l4.2.2l 

5.3. Follow-up observations of 2M1 707+64 

Follow-up observations on 2M1707H-64, also using 
the BUSCA instrument, were performed as Directors 
Discretionary Time observations on the 14th and 15th of June 
2005. The data reduction and analysis were done in the same 
way as for the 2002 and 2003 data. The relative light curves 
of 1, R and G are shown in Fig. and Fig. All three 
channels are detections according to the test. Although the 
R- and G-band data seems to consist of large random-like 
variations, the I-band light curve strongly resembles that from 
the 2002 data. Moreover, the I-band power spectrum contains 
a peak above the 0.1 % FAP power level at 3.61 + 0.2 hr which 
confirms the one we detected in the earlier data. This period 
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Figure 11. The relative light curves of the first night of follow- 
up observations of 2M 1707 +64 with BUSCA in June 2005. 1- 
band data is represented by solid lines, R-band data by dashed 
lines and G-band data by dotted lines. Error bars are shown 
only for one data point per channel. 



can also be found on night two of the G-data (P - 3.80±0.2 hr) 
at the same significance and on night one with P = 4.0 ±0.27 hr 
at a FAP of 1 %. The constancy of the period suggests that this 
is indeed a rotation period. We perhaps cannot say the rotation 
period as the object may show differential rotation, although if 
so, the spots would have to be at similar latitudes (and/or the 
differential rotation would have to be small). Due to the lack of 
a detected period in R, we could not repeat the model spectra 
procedure to investigate the source of variability. 

Further independent observations on 2M 1707 +64 were car- 
ried out at the Maidanak Observa tory, Uzbekistan, in the I - 
band in June 2005, as reported in Rockenfeller et al. (lOO&i. 
These observations were made primarily to further investigate 
the flare activity of this target, but they also strongly support the 
rotation period, detected in the Maidanak data at 3.63+0.08 hr 

5.4. Variability dependences 

One of the major goals of this survey was to study whether vari- 
ability is more common in L dwarfs than in M dwarfs. To judge 
this, we compiled a list of recently observed (1999-2005) M 
and L field dwarfs from various publications (for sources, see 
the caption of Fig.ll3>. A direct comparison of data published 
by different groups is ambiguous since different detection lim- 
its were achieved or different significance levels were used. We 
therefore only use those data which we judge to be compara- 
ble to ours. Fig.[0]shows the 1 band variability amplitudes (or 
upper limits for non-variable objects, if available) versus the 
spectral type for these samples. Fig. [H] shows variability am- 
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Figure 12. The relative light curves of the second night of 
follow-up observations on 2M1707+64 with BUSCA in June 
2005. 1-band data is represented by solid lines, R-band data by 
dashed lines and G-band data by dotted lines. Error bars are 
shown only for one data point per channel. 



plitudes and upper limits of the I-, R- and G-band data for all 
19 objects of our data set. No definite trend of the amplitudes 
with spectral type is apparent in either figure, although ampli- 
tudes of more than 0.02 mag have been found only at spectral 
types LO to L3. 

As most surveys have only been carried out in the I band, 
we use just our I band variability detections when comparing 
our results with other surveys. We then form three spectral type 
bins - early M (M0-M4, 4 objects), late M (M5-M9, 17 ob- 
jects) and L (31 objects) - and calculate the fraction of vari- 
able objects in each. We thus arrive at fractions of 0.25 + 0.25, 
0.29 + 0.13 and 0.48 + 0.12 for these bins respectively (and 
0.29±0. 12 for all M dwarfs together). Errors are derived assum- 
ing counting statistics. The difference between M and L type is 
about 1 .5 cr which corresponds to a confidence level of about 
86 %. Of particular interest is to compare only the L dwarfs 
from BJMOl with our M dwarfs because the analysis meth- 
ods and sensitivity are very similar. We now find fractions of 
0.21+0.11 for the 19 M dwarfs and 0.70 ± 0.26 for the 10 
L dwarfs, which is a 2 cr difference. Both results suggest an in- 
creased frequency of occurrence of variability among L dwarfs. 
This is supported further by the fact that the present survey is 
more likely to detect variability in the M dwarfs on account of 
the multichannel observations. 

On account of the incompleteness of further information on 
the targets as e.g. Hq. or v sin /, we cannot properly test whether 
these quantities are significantly different between variable and 
non- variable objects. 
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Figure 13. Variability amplitudes (variable objects, solid 
markers) and upper detection limits (non-variable ob- 
jects, open markers) compiled from all relevant works on 
field dwarfs in the I-band of the r e cent y ears, includ- 
ing: t h is paper. [Bailer- Jones & Mund3 (l200ll) . IClarke et all 
l|2002). Clarke et al. (2003h. lGelino et alJ l l2002h . lKoenl l l2005h . 
Martin et ah (.2001). The spectral types are slightly offset from 
the real ones in case there is more than one object of the same 
type. Amplitudes an d upper limits are always rms values. Since 
IClarkeetal.l l l2002h give peak-to-peak amplitudes for their two 
variable targets, we converted those to rms values assuming a 
sinusoidal signal (that assumption is not entirely justified, but 
diff'erences in amplitudes are less than fifteen percent). Skeletal 
markers indicate rms amplitudes of non-variable objects where 
no upper limits are available. Starred symbols stand for rms- 
amplitudes where the test gives a detection but variability is 
finally rejected. 



6. Conclusions 

We have presented multiband data of 19 M dwarfs of which 
seven show evidence for variability at a 99 % confidence level 
in at least one of four channels. We performed relative photom- 
etry along with a careful error estimation in order to achieve 
the high photometric precision needed to detect low amplitude 
variability. Amplitudes (root-mean- square of the Ught curve) 
measured from 0.0055 to 0.014 mag in I and R and from 0.014 
to 0.034 mag in G. For non-variable objects, upper limits to 
variability are estimated to lie between 0.006 and 0.05 mag. 

Using the CLEAN algorithm to form the periodograms, 
convincing periodic variability was found in LHS370 and 
2M1707+64 at periods of 5.9 ± 2.0 and 3.65 + 0.1 hours, re- 
spectively. We further claim these to be their rotation periods. 
Besides these two targets, three other objects show periodic 
variations in one channel: 2M1344-I-77 in G at 12.5 ± 2.0 hr; 
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4 6 
spectral type M 

Figure 14. Variability amplitudes (variable objects, solid markers) 
and upper detection limits (non-variable objects, open markers) of 
our data set versus M-spectral type. Channels G, R and I are repre- 
sented by triangles, squares and pentagons, respectively. Starred sym- 
bols stand for rms amplitudes where the test gives a detection but 
variability is finally rejected. 



LHS2930 in I at 13.2 ± 1 .9 hr; 2M17 14-H30 in R at 6.9 ±0.15 hr. 
But here the evidence for periodicity is less convincing. 

Various simulations were performed to determine both the 
period uncertainty and the probability of detecting an existing 
periodic signal for arbitrary time samplings. In this way we 
could confirm the high sensitivity of our data to the expected 
rotation periods of M dwarfs up to about 12 hours (although 
early M dwarfs may have longer pe riods). With the help of syn- 
thetic spectra of'A llard et aD(l200ll) the two favourite scenarios 
for variability in ultra cool dwarfs - magnetically induced spots 
and dust clouds - were studied in the cases of 2M1707H-64 and 
LHS370. It turned out that the observed variability in the for- 
mer object can be reproduced well by assuming an atmosphere 
of 2300 K, \ogg = 5.0 and a 100 K cooler spot covering 7.5 % 
of the visible star disk. Clear clouds on a dusty atmosphere or 
dusty clouds on a clear atmosphere can be ruled out. In con- 
trast to this, the origin of variability in LHS370 remains un- 
clear The lack of detectable periods for most of our targets can 
be explained by absent or rapidly evolving surface features. 

Combining the results of this and other works published in 
the recent years leads to evidence for an increased variability 
frequency in field L dwarfs when compared to field M-dwarfs. 
The fractions of variable objects are: early M type (M0-M4) 
0.25 + 0.25, late M type (M5-M9.5) 0.29 ± 0.13 and L type 
0.48 + 0.12. The number of objects in each bin is 4, 17 and 
31, respectively. The difference between M and L type is about 
1.5cr which corresponds to an 86 % confidence level. This is 
statistically not large but still hints at a more frequent variabil- 
ity in later spectral type. If we consider the homogeneous sam- 



ples of this paper and BJMOl, we find fractions of 0.21 ± 0.1 1 
(19 M dwarfs) and 0.70+0.26 (10 L dwarfs) which coiTesponds 
to a 2cr difference. 

Higher SNR data would allow us to put stronger constraints 
on variability sources, particularly if optical and infrared obser- 
vations were combined. An increased number of monitored M, 
late L and even T dwarfs could clarify the change of variability 
frequency with spectral type, and thus permit stronger conclu- 
sions. And while difficult to achieve because of the available 
telescope time, a denser and longer time sampling would make 
it possible to recover longer periods which are expected to be 
common among early M dwarfs. 

Acknowledgements. We are very grateful to Roland Gredel, the direc- 
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location of Directors Discretionary Time in 2005 and to the Calar Alto 
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Table 4. Summary of important results: A'„ is the number of reference stars, the actual and critical values (i.e. for p = 0.01), var. 

flags indicates whether a target is variable on the individual nights and in total (the latter in parenthesis) according to the test, p,ot is the 
probability of the null hypothesis, ampl. the variability amplitude (root-mean-square) or an upper limit hereto (both in magnitudes), period is 
self-explanatory (measured in hours with uncertainties obtained by the simulations, see Sect. 14.2.11 although note that theoretical values are 
larger); AR (Amplitude Ratio) is the ratio of the rms amplitude of the light curve to the noise within the target; var. states whether or not a 
channel was finally considered variable. The -I- and - flags indicate variable and non-variable objects, respectively, whereas o means that it is 
only a marginal detection (either in terms of the;^f^ test or of a consideration of all means of investigation). 



target 


SpT 


band 






var. flags 


Ptot 


ampl. 


period 


AR 


var 


2M1311-I-80 


M8 


I 


17 


38.3/54 


o (-) 


0.21 


0.012" 


— 


1.2 


- 






R 


12 


20.1/26 


-(-) 


0.065 


0.020" 


— 


1.4 


- 






G 


13 


58.5/54 


-,-(+) 


2.9 ■ 10-4 


0.088* 


— 


1.4 


- 


2M1336-H47 


M7 


I 


9 


72/95 


-,-(-) 


0.71 


0.0090" 


— 


1.3 


- 






R 


10 


66.9/95 


-,-(-) 


0.41 


0.022" 


— 


1.0 


- 






G 


9 


62.2/96 


-,-(-) 


0.61 


0.063" 


— 


0.8 


- 


2M1344-H77 


M7 


I 


33 


63.5/64 


+ (o) 


0.011 


0.0057" 


— 


2.0 


o 






R 


18 


51.8/58 


-,-(-) 


0.034 


0.016" 


— 


1.1 


- 






G 


19 


28.7/64 


-,-(-) 


0.91 


0.048" 


12.5 + 3.0? 


0.6 


- 


LHS370 


M5 


I 


14 


57.7/37 


+ 


< lO *" 


0.015* 


— 


2.4 


+ 






R 


13 


19.7/34 


- 


0.29 


0.020" 


5.9 ±2.0 


1.2 


+ 






G 


14 


16.1/37 


- 


0.65 


0.045" 


6.5 + 2.0 


1.1 


+ 


LHS2930 


M6.5 


I 


10 


80/85 


+,-(-) 


0.024 


0.026" 


13.2+ 1.9? 


2.4 


o 






R 


16 


55/85 


-,-(-) 


0.53 


0.0088" 


— 


2.1 


- 






G 


10 


56/85 


-,-(-) 


0.52 


0.063" 


— 


1.1 


- 


CTI1539-H28 


M2 


I 


29 


74.9/75 


-, (o) 


0.01 


0.0074" 


— 


2.8 


o 






R 


21 


40.9/75 


-,-(-) 


0.79 


0.0063" 


— 


0.8 


- 






G 


15 


67.0/74 


-,-(-) 


0.036 


0.014" 


— 


1.4 


- 


2M1546-H37 


M7.5 


I 


25 


58.7/75 


-,-(-) 


0.84 


0.0063" 


— 


1.2 


- 






R 


23 


105/75 


+, +, - (+) 


5 ■ 10-^ 


0.012'' 


— 


2.7 


+ 






G 


18 


75/75 


-,-(o) 


0.01 


0.025° 


— 


1.8 


- 


LHS3189 


M5.5 


I 


35 


59.3/74 


o,-(-) 


0.13 


0.0056" 


— 


1.1 


- 






R 


27 


35.6/74 


-,-(-) 


0.91 


0.0097" 


— 


0.7 


- 






G 


18 


22.2/74 


-,-(-) 


0.9995 


0.032" 


— 


0.6 


- 


2M1627-)-81 


M9 


I 


36 


56.8/78 


-,-(-) 


0.37 


0.011" 


— 


1.0 


- 






R 


31 


74.0/78 


-,-, + (-) 


0.019 


0.024" 


— 


1.4 


- 






G 


37 


41.3/78 


-,-,-(-) 


0.83 


0.051" 


— 


0.8 


- 


CTI1629-H28 


M4 


I 


28 


44.1/74 


-,-(-) 


0.64 


0.012" 


— 


1.0 


- 






R 


19 


206/74 


+, + (+) 


< 10-'' 


0.014'' 


— 


1.9 


+ 






G 


20 


86.8/73 


+,-(+) 


4 ■ 10-4 


0.022* 


— 


1.5 


- 


2M1707-t-64 


M9 


I 


37 


143/54 


+, + (+) 


< lO-*" 


0.012'^ 


3.65 ±0.1 


2.5 


+ 






R 


32 


62/54 


+, + (+) 


0.001 


0.014* 


3.7 + 0.1 


1.7 


+ 






G 


25 


50/54 


+ (-) 


0.022 


0.029" 


3.3 + 0.1 


1.8 


+ 


CTI1709-H27 


M5.5 


I 


43 


120/80 


+ (+) 


< lO-*" 


0.0062* 




2.1 


+ 






R 


26 


277/80 


+, + (+) 


< lO-*" 


0.011'' 




2.2 


+ 






G 


29 


140/78 


+, +, + (+) 


< lO-*" 


0.014* 




1.6 


+ 


2M1714-H30 


M6.5 


I 


32 


30.0/77 


-,-(-) 


0.99 


0.0080" 




0.8 








R 


31 


107/77 


+, + (+) 


< lO-*" 


0.012* 


6.9 + 0.15 ? 


1.9 


+ 






G 


27 


144/77 


+, + (+) 


5 • lO-*" 


0.034* 




1.8 




LHS3307 


M5 


I 


49 


67.5/84 


-,-(-) 


0.13 


0.052" 




> 3 








R 


50 


104/84 


+,-(+) 


0.0001 


0.012* 




1.2 









G 


47 


57.7/84 


-,-(-) 


0.41 


0.031" 




0.8 




2M17504-44 


M7.5 


I 


31 


84.3/92 


+ (-) 


0.038 


0.0075" 




1.1 








R 


36 


62.4/88 


-,-(-) 


0.36 


0.018" 




1.3 








G 


33 


69.7/80 


o, - (-) 


0.061 


0.036" 




1.2 




LHS3339 


M6 


I 


32 


116/75 


-,-(+) 


< lO-*" 


0.0061* 




1.7 


+ 






R 


31 


67.6/74 


-,-(-) 


0.033 


0.0083" 




1.5 








G 


20 


119/75 


-,-(+) 


< lO-*" 


0.021* 




1.3 


+ 
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target 


SpT 


band 




X Ixq 


var. flags 


Ptot 


ampl. 


period 


AR var. 


2M1757+70 


M7.5 


I 


50 


48.6/60 




0.097 


0.0061" 




4.7 - 






R 


39 


54.8/60 




0.03 


0.0090" 




1.5 - 






G 


29 


50.7/60 


+ (-) 


0.066 


0.025" 




1.8 - 


CTI1801+28 


M2.5 


I 


31 


49.0/84 




0.73 


0.0061" 




1.1 - 






R 


31 


73.7/81 




0.039 


0.0055" 




1.2 - 






G 


27 


63.2/81 




0.18 


0.015" 




1.1 - 


LHS3376 


M4 


I 


44 


135/50 


+ 


< 10-6 


0.052* 




> 3 - 






R 


25 


27/50 




0.43 


0.029" 




1.3 - 






G 


40 


45/50 




0.029 


0.050" 




1.9 - 



amplilLidc column: " upper-limil-ims; '' lighl curve mis 



